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Abstract 



Fits of pion-nucleus potentials to large sets of pionic atom data reveal 
departures of parameter values from the corresponding free ttN parameters. 
These medium effects can be quantitatively reproduced by a chiral-motivated 
model where the pion decay constant is modified in the medium or by in- 
cluding the empirical on-shell energy dependence of the amplitudes. No con- 
sistency is obtained between pionic atoms and the free ttN interaction when 
an extreme off-shell chiral model is used. The role of the size of data sets is 
briefly discussed. 
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I. INTRODUCTION 



Renewed interest in pionic atoms in general, and in the s-wave part of the pion-nucleus 
potential in particular, stems from three recent developments. The first is the experimental 
observation of 'deeply bound' pionic atom states in the (d, 3 He) reaction [1-4], the existence 
of which was predicted a decade earlier [5-7] . The second is the very accurate measurements 
of the shift and width of the Is level in pionic hydrogen [8] and in pionic deuterium [9,10] 
which leads to precise values for the s-wave nN scattering lengths (see also Ref. [11]). The 
third development is the attempt to explain the 'anomalous' s-wave repulsion [12] in terms 
of a density dependence of the pion decay constant [13], or very recently by constructing the 
Ti N amplitude near threshold within a systematic chiral perturbation expansion [14] and in 
particular imposing on it gauge invariance [15,16]. 

The so-called anomalous repulsion of the s-wave pionic atom potential is the empirical 
finding, from fits of optical potentials to pionic atom data, that the strength of the repulsive 
s-wave potential inside nuclei is typically twice as large as is expected on the basis of the 
free ttN interaction. This enhancement is, to a large extent, due to the value of the potential 
parameter b\, the in-medium isovector s-wave ttN amplitude. The possible modification in 
the medium of this amplitude is the topic of the present work. One may define medium effects 
in this context as unexpected results obtained from analyses of experimental data, usually 
with rather simple models, based in the present case on the free pion-nucleon interaction. 

Section II presents briefly the standard pion-nucleus potential [17] where the various 
parameters are defined. Section III shows results of fits of the standard potential to a large- 
scale ('global') set of data where medium effects are presented in a quantitative way. Section 
IV shows similar results for the model due to Weise [13] where the pion decay constant is 
assumed to be modified in the medium. An alternative approach where the medium effects 
are essentially due to the energy dependence of the amplitudes is discussed in section V 
both for a chiral (off shell) approach and for an empirical (on-shell) approach. In section 
VI the statistical significance of the results and their relation to the extent of the data base 
are discussed together with consequences of constraining parameter values. Section VII is a 
summary. 



The optical potential V opt used in the Klein-Gordon equation is of the standard form due 
to Ericson and Ericson [17], 



II. THE PION-NUCLEUS POTENTIAL 
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where 



ai (r) = 4tt(1 + A)-i{ Co [ Pn ( r ) + pp ( r )] + Cl[pn ( r ) _ pp ( r )]} ; ( 4 ) 

a 2 (r) = 4tt(1 + ^)- 1 4C p„(r) Pp (r). (5) 

In these expressions p n and p p are the neutron and proton density distributions normalized 
to the number of neutrons N and number of protons Z, respectively, /i is the pion-nucleus 
reduced mass and M is the mass of the nucleon, q(r) is referred to as the s-wave potential 
term and a(r) is referred to as the p-wave potential term. The latter will not be discussed 
here further as there are essentially no non-trivial medium effects with this term. The 
function bo (r) in Eq. (2) is given in terms of the local Fermi momentum kp(r) corresponding 
to the isoscalar nucleon density distribution: 

b (r) = b -^-(bl + 2b 2 1 )k F (r), (6) 

where bo and b\ are minus the pion-nucleon isoscalar and isovector effective scattering 
lengths, respectively. The quadratic terms in b and b\ represent double-scattering mod- 
ifications of bo- In particular, the b\ term represents a sizable correction to the nearly 
vanishing linear b term. 

The nuclear density distributions p p and p n are an essential part of the pion-nucleus 
potential. Density distributions for protons can be obtained by unfolding the finite size of 
the charge of the proton from charge distributions obtained from experiments with electrons 
or muons. This leads to reliable proton densities, particularly in the surface region, which 
is the relevant region for producing strong-interaction effects in pionic atoms. The neutron 
distributions are, however, generally not known to sufficient accuracy and we have therefore 
adopted a semi-phenomenological approach that covers a broad range of possible neutron 
density distributions. The feature of neutron density distributions which is most effective 
in determining parameter values of the potential is the radial extent, as represented, for 
example, by r n , the rms radius of the neutron density. Other features such as the detailed 
shape of the distribution have only minor effect on values of potential parameters, although 
they do affect the quality of fits. For that reason we chose the rms radius as the prime 
parameter in the present study, and focus attention on values of r n — r p , the difference 
between the rms radii. A simple parameterization adopted here that is easy to relate to the 
results of relativistic mean field (RMF) calculations [18,19] is 

N ~ Z 

r n -r p = a - + rj. (7) 

In order to allow for possible differences in the shape of the neutron distribution, the 'skin' 
and the 'halo' forms of Ref. [19] have been used. Note that the results of RMF calculations 
are described very well [19] by using a— 1.5 fm and rj— —0.03 fm. This value of a is probably 
an upper limit as other sources of information suggest values close to 1-1.2 fm. 
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III. RESULTS FOR THE STANDARD POTENTIAL 



Results obtained from the use of the standard potential of the previous section, here 
marked as 'conventional', with 100 data points along the periodic table, from 20 Ne to 238 U, 
are shown in Fig. 1. These data include also the deeply bound Is states in 205 Pb [3] and 
in 115 > 119 . 123 g n [4]. The figure shows that for a— 1.5 and for the 'skin' type of density a very 
good fit to the data is obtained (x 2 per point of 1.7) but the resulting value of b\ is more 
repulsive than the free ttN value by at least two error bars, or three error bars if the more 
reasonable value of a— 1.25 is adopted as representing the average behaviour of rms radii 
of neutron distributions. This enhanced repulsion is a clear indication of medium effects. 
Other medium effects in the potential are likely to be present in the phenomenological two- 
nucleon absorption terms Bq and Co which also have real dispersive parts. These are more 
difficult to handle and we will not discuss these here except for noting that the real parts of 
B and C are expected not to exceed in absolute values the corresponding imaginary parts, 
respectively. The anomalous s-wave repulsion observed with the conventional potential is 
due to the combined action of the too repulsive b\ (which enters also quadratically, see Eq. 
(6)) and of Rei?o which is also found to be too repulsive compared to the above expectations 
(see also section VI below). 



IV. MEDIUM-MODIFIED PION DECAY CONSTANT 

The in-medium s-wave interaction of pions have been discussed recently by Weise [13] 
in terms of partial restoration of chiral symmetry in dense matter. Since b\ in free-space is 
well approximated in lowest chiral-expansion order by the Tomozawa- Weinberg expression 
[20] 

* = = -<>■«< W 

then it can be argued that b\ will be modified in pionic atoms if the pion decay constant f n 
is modified in the medium. The square of this decay constant is given, in leading order, as 
a linear function of the nuclear density, 

f: 2 = fl - — 2 p (9) 

mi 

with a the pion-nucleon sigma term. This leads to a density-dependent isovector amplitude 
such that bi becomes 

bi(p) = (10) 
Kyj l-2.3p V ; 

for cr=50 MeV [21] and with p in units of fm -3 . This model was found [22,23] to be very 
successful when tested against large data sets. 

Figure 2 shows results similar to those of Fig. 1 but with bi as given by Eq. (10). 
Only the results for the 'skin' shape of the neutron densities are plotted as the 'halo' shape 
always produces poorer fits to the data. It is seen from the figure that the quality of fits to 
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the data is the same as for the conventional model, but the values of b\ have shifted now 
such that for a ~ 1.25 they agree with the free 71N value, thus indicating that the density 
dependence of the decay constant could account for the medium modification of b±. Using 
the present model the parameter Re_B is found to be considerably less repulsive than in the 
conventional model, and within errors in agreement with expectations [22]. 

V. ENERGY-DEPENDENT AMPLITUDES 

Recently, Kolomeitsev, Kaiser and Weise [15] have suggested that pionic atom data could 
be reproduced using a pion optical potential underlain by chirally expanded 7riV amplitudes, 
retaining the energy dependence of the amplitudes bo(E) and b\(E) for zero- momentum 
(q=0) pions in nuclear matter in order to impose the minimal substitution requirement 
E — > E — V c , where V c is the Coulomb potential. This has the advantage of enabling one 
to use a systematic chiral expansion as an input [14], rather than singling out the leading 
order term Eq. (8) for b\. 

The chiral expansion of the nN amplitudes for q = at the two-loop level is well 
approximated by the following expressions [14,15]: 

*( 1 + £)w^(^ + Sg), 

qa is the nucleon axial-vector coupling constant, qa = 1.27, (3 and 7 are tuned to reproduce 
the threshold values &o( m 7r) ~ and bi(m n ) = — 0.0885to;oo2i [8] respectively. For b , 

in view of the accidental cancellations that lead to its near vanishing we limit our discussion 
to the f~ 2 term in Eq. (11), therefore choosing /3 = am~ 2 . 

Implementing the minimal substitution requirement in the calculation of pionic atom 
observables, the constant parameters 60,1 °f the conventionally energy independent optical 
potential have been replaced in our calculation [24] by 

Mr) = 60,1 " So^ReB + V c (r)) , (13) 

where 80,1 = dbo^(E)/dE is the appropriate slope parameter at threshold, KeB is the (real) 
binding energy of the corresponding pionic atom state and V c (r) is the Coulomb potential. 
The constant fit parameters fe ,i are then expected to agree with the corresponding free nN 
threshold amplitudes if the energy dependence is indeed responsible for the renormalized 
values found in conventional analyses. 

In addition to the above 'chiral' energy dependence for off-shell q = pions we also 
present results for the empirically known on-shell nN amplitudes, when the pion energy 
E and its three- momentum q are related by E 2 = m 2 + q 2 . This choice corresponds to 
the original suggestion by Ericson and Tauscher [25] to consider the effect of energy depen- 
dence in pionic atoms. Ericson subsequently [26] pointed out that, for strongly repulsive 
short-range iViV correlations, the on-shell requirement follows naturally from the Agassi-Gal 
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theorem [27] for scattering off non-overlapping nucleons. The corresponding itN amplitudes 
are denoted as 'empirical' and are taken from the SAID data base [28]. 

Figure 3 shows results for the two kinds of energy-dependence mentioned above. The left 
panels show that implementing the chiral off-shell energy dependence leads to very significant 
deterioration in the fit to the data both for pionic atoms (lower part) and for the free nN 
interaction (upper part). Note also that the minimum in the x 2 curve occurs at unacceptably 
large value of the parameter a. In contrast, the right hand side of the figure shows good 
agreement both for pionic atoms and for the free nN interaction, in accordance with the 
original suggestion made by Ericson and Tauscher [25]. Perhaps it is not too surprising 
that the off-shell q=0 nuclear matter approach is inadequate when the pion-nuclear optical 
potential that generates pionic-atom wavefunctions is being considered. 

VI. STATISTICAL CONSIDERATIONS AND DATA SETS 

The results discussed so far have all been based on the 'global 3' data set of Ref. [19] 
which consists of 100 data points. Any conclusion regarding medium effects must obviously 
be linked to the uncertainties in the extracted parameter values, and these invariably depend 
on the size of the data sets used. In addition the uncertainties may depend critically on 
assumptions made in the analysis such as assigning fixed values to some parameters. These 
points are demonstrated in Table I. 

The different columns of the table refer to various data sets, ranging from 120 points for 
nuclei between 12 C and 238 U to just 20 points where the deeply bound states provide half of 
the data. The top half of the Table is for unrestricted x 2 fits whereas the lower part is for 
corresponding fits with the parameter KeB set to zero. All the results in this table are for 
a conventional potential, i.e. without explicit medium effects such as density-dependence 
or energy-dependence of the amplitudes, and for neutron densities calculated for a— 1.5 
(Eq. (7)). From the top half of the table it is seen that the values of b are essentially in 
agreement with the free nN value, within their rather large uncertainties and that the values 
of b\ disagree with the corresponding free nN value only for the two large data sets since the 
uncertainties become too large for the two smaller data sets to reach such a conclusion. This 
distinction becomes even sharper when the uncertainty due to the neutron distributions is 
also included [19]. Another feature of the top half of the table is that the parameter ReB 
is distinctly different from zero, although its uncertainty is not small, only when the large 
data bases are being used. For the two smaller data sets its value is consistent with zero. 

Some recent papers insist on using the value of Re_B =0 [2-4] and the lower part of the 
table addresses the consequences of such a constraint. The first consequence is that the 
uncertainty of the parameter 6 becomes much smaller than before and its values, in all 
cases, deviate sharply from the free nN value. This conclusion remains valid also for all the 
types of medium effects for b\ considered in the present work. The second consequence is 
that for the large data sets the constraint Rei?o=0 indeed leads to a significant deterioration 
in the quality of fits, as evidenced by the increase of x 2 ■ Recall that the natural unit for 
such increases is the value of x" 2 '/F, the x 2 P er degree of freedom. For the large data sets 
the increase is 10 such units, which is very significant. Note also that the values of b\ do 
not depend on whether or not R,eB is set to zero. Similar results are obtained for the other 
models for b\ discussed above, but with values of b\ in agreement with what is displayed 
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in the figures. It is concluded that analyses of reduced data sets with the constraint of 
Rei?o=0 must lead to unreliable results. In fact, analysing only the deeply bound states 
for 115 ' 119 . 123 g n anc l 205 Pb, we obtain excellent fits to the data with values of b\ anywhere 
between —0.07 and —0.13 m^ 1 . 

VII. SUMMARY 

We have shown that global fits to large sets of data on strong interaction effects in pionic 
atoms in terms of a theoretically-motivated optical potential lead to very good description of 
the data (x 2 /N ~1.7) with well-determined s-wave isovector amplitude that differs from the 
corresponding free nN amplitude by three standard deviations. This difference indicates 
modification of the interaction in the medium and it was shown that a chiral-motivated 
model where the pion decay constant is modified in the medium is capable of reproducing 
the medium effects. An alternative approach where the empirical on-shell energy dependence 
of the 7rN amplitude is included within the minimal substitution E — > E — V c was also shown 
to be successful. In contrast, the fully off-shell chiral energy dependence of the amplitudes 
fails badly. The importance of unrestricted and unbiased fits to large scale data sets has 
been demonstrated. 

I wish to thank A. Gal for many fruitful discussions throughout the whole project. This 
research was partially supported by the Israel Science Foundation grant No. 131/01. 
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FIGURES 




FIG. 1. Values of x 2 f° r 100 data points (lower) and best-fit values of 61 (upper) as function of 
the neutron density parameter a. 'exp' marks the experimental value for the free irN interaction 
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FIG. 2. Values of \ 2 f° r 100 data points (lower) and best-fit values of 61 (upper) as function of 
the neutron density parameter a. 'exp' marks the experimental value for the free ttN interaction 
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FIG. 3. Values of x 2 for 100 data points (lower) and best-fit values of b\ (upper) as function of 
the neutron density parameter a. 'exp' marks the experimental value for the free 7riV interaction. 
Left panels, for the chiral energy dependence of &o and b±; right panels, for the empirical energy 
dependence. 
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TABLES 



TABLE I. Dependence of parameter values and uncertainties on the extent of data base and 
on assumptions made regarding Rei?o- Deterioration in the quality of fits is measured by p^p, the 
increase of \ 2 in units of % 2 per degree of freedom. 



Values of r n 


are for a=1.5 








data 
points 


'global 2' 

12 C to 238 U 

120 


'global 3' 
20 Ne to 238 U 

100 


light N = Z 
+ light N > Z 
Is only 
22 


light N = Z 
+'deep' 
Is only 
20 


x 2 

X 2 /F 

him- 1 ) 

ReB (m~ 4 ) 

ImB (m~ 4 ) 


237 
2.1 
0.000(6) 
-0.101(3) 
-0.085(30) 
0.049(2) 


171 
1.8 
-0.001(7) 
-0.098(3) 
-0.082(30) 
0.052(2) 


54 
3.0 
-0.009(17) 
-0.095(13) 
-0.048(72) 
0.049(2) 


35 
2.2 
-0.016(13) 
-0.094(7) 
-0.017(60) 
0.051(2) 


x 2 

X 2 /F 

Ax 2 
X 2 /F 

him- 1 ) 

ReB (m- 4 ) 

ImBo(m~ A ) 


259 

2.3 

9.6 
-0.018(1) 
-0.102(3) 
0. (fixed) 
0.048(2) 


190 

2.0 

10.6 
-0.019(1) 
-0.099(3) 
0. (fixed) 
0.051(2) 


56 
2.9 
0.7 
-0.020(3) 
-0.093(12) 
0. (fixed) 
0.048(2) 


35 
2.1 


-0.020(2) 
-0.094(7) 
0. (fixed) 
0.050(2) 



'deep' refers to deeply bound Is states in 115 > 119 > 123 Sn and 205 Pb. 

b = -0.00011^021 rn' 1 and b x = -0.08851^02? rn' 1 for the free vriV interaction. 
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